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ABSTRACT 

The linear boundary layer solution (Chang, 1973b), 
which includes the temporal acceleration, is utilized 
for a CISK parameterization to study tropical lower tro- 
pospheric waves. A one-level equation is derived for 
the top of the boundary layer and is solved numerically 
for its eigenvalues. Calculations are made for varying 
boundary “lay er depths and vertical heating gradients. 
Growth rates, periods and eigenfunctions for the lowest 
meridional modes of synoptic-scale disturbances are 
analyzed. Results indicate that CISK heating does not 
prefer any finite wavelength. 
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INTRODUCTION 



I . 

Observational evidence of tropical lower tropospheric 
waves has been extensively collected and studied during 
the past several years. Although there is still some doubt 
as to the detailed characteristics of these disturbances, 
it is generally observed that these waves propagate west- 
ward, exhibit a period of about four to five days in their 
meridional wind component, and may fall in one of the two 
wavelength ranges: around 2000-4000 km, or approximately 

8000 km. As a direct result of the observations, meteoro- 
logists have attempted to explain theoretically the 
mechanism and structure of these tropical disturbances. 

Three different hypotheses have thus far been postu- 
lated: 1) lateral forcing by mid-latitude systems, 2) 

barotropic instability, and 3) thermal forcing. A concise 
but excellent description of these theories is found in 
Yanai (1971). To summarize his discussion briefly, various 
investigators have found that all three hypotheses may be 
responsible for tropical wave growth. 

Widespread attention has been given the thermal forcing 
mechanism during the past decade. Ever since the inception 
of the so-called CISK (Conditional Instability of the 
Second Kind) theory, first proposed by Charney and Eliassen 
(1964) and Ooyama (1964) to explain hurricane growth, many 
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investigators have sought to extend its applicability to 
the larger scale tropical waves. Basically their CISK 
models incorporate thermal forcing by parameterizing the 
latent heat release from small scale convection in terms 
of the boundary layer convergence of the large scale flows. 
In addition, Charney and Eliassen’s geos trophic/Ekman 
boundary layer formulation (geostrophic vorticity and ver- 
tical velocity are proportional at the top of the boundary 
layer) was commonly used. 

The primary approach in many of the previous studies 
was to perform linear stability analysis of a system of 
equations in order to determine a preferred wave scale - 
that scale which has the largest growth rate when other 
conditions remain constant. Most of these earlier linear 
studies used either a geo s tr ophi c/Ekman boundary layer or 
a surface drag representation for the lowest level in the 
model. ' Although the linearization process was done to 
simplify an analysis, a question generally arose in the 
Ekman-type boundary layer models near the equator. In this 
region the Coriolis force in the equations of motion 
approaches zero and therefore the geostrophic relation may 
no longer be valid. It is pointed out by Holton et al 
(1971) that a transition zone exists near the equator where 
the Ekman boundary layer structure of the more poleward 
latitudes changes rather drastically into a Stokes-type 
boundary layer in the equatorial region. Generally, this 
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transition zone has not been accounted for by most of the 
investigators . 

For some of the earlier CISK wave models, a major dis- 
crepancy exists. Yamasaki (1971) and Murakami (1972) 
found that preferred wave scales on the order of several 
thousand kilometers exist in their models. However, in 
contrast to these results, Chang (1971) and Williams and 
Robertson (1973) showed that a distinct maximum growth rate 
for the observed tropical waves did not occur. 

Because of this discrepancy, it remained inconclusive 
as to the effect CISK-type heating had upon tropical waves 
in the lower troposphere. The purpose of this research is 
to use an improved boundary layer representation within the 
scope of linear theory to study the effect of CISK upon 
tropical wave growth. The linear boundary layer solution 
(Chang, 1973b) which includes the temporal acceleration is 
used for the CISK parameterization. By neglecting the small 
temperature fluctuations in the thermodynamic equation, 
which is generally- valid for weak synoptic- scale disturbances 
in the lower tropical atmosphere, a one-level equation is 
derived for the top of the boundary layer. This equation is 
similar to the free atmosphere equation in Chang (1973a). 

An eigenvalue problem is then formulated and the growth rates 
of the lowest meridional modes are obtained numerically. 
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II. THE MODEL 



A. BASIC EQUATIONS 

The basic equations of the model are the linearized 
Boussinesq primitive equations for an equatorial beta- 
plane: 



9u* 

3t* 


3p* . „ 3 u * 

* - 3x* + By v + K . ,2 ' 

dz* 


(1) 


3 v * 
3t* 


3p* 0 * „ 3 2 v* 

' - 3y* - 6y *“* + K \ ,2 ' 

J 3z * 


(2) 


3 u * 
3x * 


3v* 3w* _ 

+ 3y* + 3z* 0 ' 


(3) 



where 

x*,y*,z*,t* 



zonal, meridional, vertical and time 
coordinates , respectively 



u* (x* , y * , z* , t* ) = 
v* (x* , y * , z* , t * ) = 
w*(x*,y*,z*,t*)= 
P* (x* ,y * / 1 * ) = 
K* « 



zonal velocity 

meridional velocity 

vertical velocity 

pressure divided by density 

vertical eddy viscosity coefficient, 
assumed constant 



$ = local change of Coriolis parameter, F, 

with respect to latitude, y, at the 
equator; hence, equal to 20,/a with 0, 
being the angular velocity of the earth 
and a being the radius of the earth. 



In addition, the thermodynamic energy equation 



9t* 

3t* 



u 1 



8t* 
9x * 



+ 



3y * 



w* S * 



Q* , 



(4) 
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where 



temperature 



T* = 



S* = 



q dT* 

static stability parameter ■ 

where g is the acceleration c p 
due to gravity and c is the 
specific heat of dry^air at constant 
pressure ) 



= diabatic heating rate 



is used . 

The local change and horizontal advection terms in (4) 
can be neglected considering the magnitude of each term 
(Wallace, 1971; Holton, 1972). The thermodynamic equation 
then becomes 



w *S * = Q* 



(5) 



This equation implies that the diabatic heating is balanced 

by the adiabatic vertical motion field. 

The basic equations are non-dimensionalized using a 

1/2 

gravity wave speed, c = (gH) , where H is the scale height 
and the depth scale, and the Rossby parameter, 3- These two 

parameters are used to derive the non-dimensional units: 

- 1/2 1/2 
1) time, T = ( cB) , 2) horizontal length, L = (c/3) , 

2 

3) horizontal velocity, c = L/t , 4) pressure, 7T = c /L , 

2 

5) temperature, 0 = c /R, and 6) diabatic heating rate, 

3 , 

<}> - c c / RL . 

P 

The non-dimensional equations are: 



3u 

at 




yv + K 




( 6 ) 
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90 

3t 



yu + K 



(7) 



3P 

3y 



2 /v 

9 v 



9z ‘ 



A rv A r\ A 

9u d v dw 

9x + 9y + 9z 



0 , 



( 8 ) 



/\ /v 

w S - Q 



(9) 



where the non-dimensional static stability parameter is 

RT dT 

s = + — . 

c dz 
P 



are assumed to be propagating waves 
X, and frequency V: 



Solutions for (6) - (9) 

with zonal wave number 




Substituting (10) into (6) - 



i Vu = 


- i Xp 


+ yv + K 


3 u 
9z 2 


i V v - 


9p 

9y 


- yu + K 


9 2 v 
9z 2 


iXu + 


+ 

9y 


iiL = o 
9 z ° ' 




w S = 


Q • 







u (y , z ) 
v (y , z ) 
w (y , z ) 

P (Y ) 

Q (y / z) 

(9) , the equations become 

(ID 

( 12 ) 

(13) 

(14) 



e iUx+Vt) .(10) 
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By eliminating u, v and w in (11) - (14), a single 

differential equation in p results: 



9y 2 



- [- 



2y 



Y V 
2 



+ [1 

9y l v 



, 2 2 . 
(y + v ) 

,2 2 „ 

(y - v ) 



X 2 ] p 



i (Y 



v 2 ) 



V 



_1 (Q) 
9z V 



(15) 



In this equation the only vertically differentiated term 
is the heating term on the right-hand side. This implies 
that the system of equations is vertically coup] ed only 
through diabatic heating. Hence, if the vertical profile 
of Q is specified, (15) can be applied as an ordinary dif- 
ferential equation with y as the only independent variable 
at any level above the boundary layer. 

B. THE CISK EQUATION 

The CISK process is parameterized so that the diabatic 
heating function, Q(z), is related to the vertical velocity, 
w^, at the top of the boundary layer (z^) ; 

Q ( z ) = H (z ) * w T , (16) 

where r\ ( z ) is a proportionality function which specifies 
the vertical profile of heating. Notice that (16) also 
implies negative heating when downward motion occurs at 
the top of the boundary layer. This is difficult to avoid 
due to the separation of the x-dependence of the solutions. 

The vertical motion at the top of the boundary layer 
is expressed by 
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( 17 ) 



where 



8 p 9p . 

't = S 1 ^7 * 8 2 3 ; * AR 3 P 



R, = B C - B C 
1 + + - - 



R_ = 



2 2Ky + + 



( B D + B D ) 



+ B + C + [(y + V) 1 - 2 iyz T B + A + 2 ] 



- 1 2 
+ B_C_[(y-V) + 2iaz T B A_ ] 



R_ = 



3 2Ky + + 



(B D - B D ) 



and 



+ B + C + [(y+V) 1 - 2 ±yz t B + a + 2 + A) 



+ B C_[(y-V) ^ - 2iaz T B A_ 2 - A] 



A 



exp (- yz t ) 



A_ = exp(-az T ) 

B + = | ( y +v ) _1 <! " a + 2)_1 

B = f(y-V)" 1 (1 - A_ 2 ) _ 1 

1 2A + " 2 1 

C + Z T Y + Y A+ ^ Z T + y ^ 

1 2 1 

C = Z m - - + - A (z + -) 

T a a T a 

D = - (2A -1-A 2 ) + 2z A [— ( 1 - A )- 

Y ' 



z A ] 
T + 



D = — (2A -1-A 2 ) + 2 z A [-(1-A )- z A ] 

- a 2 - T-a - T - 
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and 



2 i (v + y) 

1 K 

2 _ i(V-y) 

K 

Equation (17) is derived by integrating the total boundary 
layer convergence, using the solutions in Chang (1973b). 

This boundary layer solution includes the temporal accelera- 
tion and, therefore, is frequency dependent. If the bound- 
ary layer depth becomes infinite (z -> 00 ) , a singularity in 
exists at the critical latitude where y = V (Holton et al, 
1971). If is finite, the solution may possess quasi- 

singular behavior around this latitude depending upon 
and the mode of symmetry about the equator (Chang, 1973b). 

The heating profile, q ( z ) , is specified such that 3q/3z. 
varies from IS to 2S at the top of the boundary layer. Con- 
sequently (15) can then be applied at z^. [Notice that 
(13) and (16) require that q = S at the top of the boundary 
layer. For a lapse rate of 6°C km ^ and w =1 cm sec ^ the 
magnitude of q in turn produces a heating rate of approxi- 
mately 3°C day Because the slope of the heating func- 

tion is defined, it is implied that the heating function 
varies smoothly in the vicinity of z = z^. Furthermore, a 
positive slope implies a heat source above the top of the 
boundary layer. Finally, to reconcile the previous 
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assumption that the pressure perturbation is independent 

of height for z < z , it is assumed that the smaller but 

— T 

finite heating in the uppermost layer of the boundary 
layer acts in such a way that p remains constant in height. 

A second-order homogeneous differential equation in p 
is obtained by substituting (16) and (17) into (15): 



E 




F 




G p 



0 , 



(18) 



where E, F, and G are functions of y and V: 



E = 1 + 



. , 2 

i (y - v ) 



vs 



dp 

dz ' 



F = - 



2 y 



(y 2 - v 2 ) 



. , 2 2 , 
1 (y - v ; 



vs 



R , 



dp 

dz ' 



r - ^ (y 2+ v 2 ) ,2 

G = V ,2 2 " A + 

(y - v ) 



. , 2 2 . 

i (y - v ) 

vs 



Xr. 



dp 

dz 
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III. SOLUTION PROCEDURE 



To solve (18) numerically, a central difference scheme 

\ 

in terms of the independent variable, y, is utilized. The 
difference expression is 



A channel-flow boundary condition is assumed such that the 
perturbations vanish at the boundaries; that is, 



The difference equation, (19), is solved for its eigen- 
value, V, using an iterative procedure as described by 
Stone (1970). First, for a given value of X, z / K and 
3q/9z , an initial guess for V is made. Then by satisfying 
one of the boundary conditions in (20), the pressure per- 
turbation at the first grid point is set to zero; that is, 
p = 0. Next, the amplitude of the eigenfunction at the 
second grid point is arbitrarily chosen as unity. Using 




0 (19) 



i = 2, 3, 4 



, . . . , 



j-i, 



where J is the total number of grid points and 



Ay 



width between the two boundaries 
J-l 




P 



J 



0 



( 20 ) 
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these first two values for the eigenfunctions, the numeri- 



cal solution process successively solves for the remaining 
values of the eigenfunctions: , p^,...Pj. 

In general, the guess for V is incorrect and the calcu- 
lated value for p is therefore non-zero, thereby violating 

u 

the second boundary condition of (20). Consequently, the 
calculation of p is repeated, using a value of V which dif- 

tJ 

fers from the initial guess by 1%. The two values of p are 

vJ 

then used to perform a linear extrapolation from the two 
values of V to a new value of V which would make p zero if 
it depends on V in a linear fashion. The new value of V is 
then used to reiterate the whole procedure, and iterations 
are repeated until a value of V accurate to 0.01% is ob- 
tained. In summary the process is as follows: 

N+l N-l Av N - 1 , . 

v ' v - S? ' > 1211 

U 



where N is the iteration step. 



Av 

Ap j 



N N-l 

V - V 



N N-l 

P J " P J 



The upper limit for N is prespecified as 25 in almost all 
the computer runs. Convergence to a solution generally 
occurs well before this upper bound but, if convergence is 
not achieved by N=25, that particular run is terminated 
and a new set of input values is automatically provided to 
maximize computer efficiency. 
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Two external parameters are varied in the following 
range: z = 1 km to 4 km and 3r|/3z = is to 2S. The solu- 

tions are obtained for both the symmetric (pressure is 
symmetric about the equator) and the asymmetric (pressure 
is asymmetric about the equator) cases. 
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IV. RESULTS 



A. FREE MODE 

Solutions for the "free" mode, when the right-hand 
side of (15) is set to zero, are obtained as a test for 
the model. Zero growth rates (infinite e-folding times) 
result for all wavelengths. 

In Figure 1, the results for wave period are shown. 
Notice that wavelengths of the synoptic scale range (2000 
-10,000 km) have periods between approximately three and 
seven days for the asymmetric case and five and ten days 
for the symmetric case. In general, the longer waves 
have higher frequencies than the shorter waves in both 
cases. These free mode frequencies correspond quite 
closely to those described by the barotropic Rossby wave 
frequency formula 



v = x/ ( A 2 + a 2 ) , 

where £ = Tf/D, D being the width of the channel, and the 
approximated Matsuno (1966) wave-frequency formula 

V = k/(k 2 + 2n + 1) . 

Here k is the non-dimensional zonal wave number and n is 
the meridional mode number. 
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B. GROWTH RATE 



Figures 2 through 5 indicate the growth rates for both 
symmetry cases. The important result shown by all four of 
these figures is that the synoptic scale waves have essen- 
tially the same growth rate. Thus, there exists no pre- 
ferred wavelength in the synoptic scale range - a result 
which agrees with those of Chang (1971) and Robertson (1973). 

For each set of external parameters, two different solu- 
tions are generally found: a long-wave (LW) mode, one that 

possesses greater growth rates in the long- wave portion of 
the synoptic scale; and a short-wave (SW) mode, one that 
favors short-wave growth. Both modes of the symmetric case 
generally show larger growth rates (faster e-folding times) 
than those of the asymmetric case, although for shorter 
wavelengths this difference rapidly decreases. 

An increase in the slope of the heating function, r| , 
increases the growth rates of the waves. For a boundary 
layer depth of one kilometer, both modes of the asymmetric 
and symmetric cases indicate an approximate doubling of 
their growth rates with an increase of dr\/dz from Is to 2s 
(compare Figures 2 and 4) . For = 2 km, the same increase 

in the vertical gradient of heating produces an even larger 
change in the growth rates - generally a 125% increase 
results (see Figures 3 and 5) . 

Likewise, a deeper boundary layer produces larger growth 
rates. This effect is especially noticeable for the LW 
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mode of the asymmetric case. Growth rates for this mode 



approximately triple when the boundary layer deepens from 
one to two kilometers, but increase more slowly with 
deeper boundary layers (results not shown) . The SW mode 
for the asymmetric case, as well as the two symmetric 
case modes, exhibits similar but less pronounced increases 
in growth rates when the boundary “layer depth increases 
from one to two kilometers. 

A rather interesting result is that the symmetric SW 
mode for shallow boun da r y -1 ay er depths 1.3 km for a 

wavelength of 4500 km), a low vertical heating gradient 
(3r|/3z = IS) and longer (>^4000 km) wavelengths has larger 
growth rates than the LW mode (Figure 4) . But for deeper 
boundary layers the opposite occurs; that is, the growth 
rate of the symmetric LW mode is greater than those of the 
SW mode (Figure 5) . Thus, a simple relationship between 
the growth rates of the symmetric modes and the boundary- 
layer depth is not apparent. On the other hand, Figures 
2 and 3 indicate the asymmetric SW mode has larger growth 
rates than those of the LW mode for both bounda ry- layer 
depths . 

C. PERIOD 

As was mentioned earlier the periods of the free modes 
(Figure 1) were determined as a test case. Figures 6 
through 9 illustrate the effects of heating and boundary- 
layer depth upon wave period. Note that the periods for 
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the LW modes of both the asymmetric case and symmetric 
cases very closely resemble those of the corresponding 
free modes. This result is expected, because the LW modes 
generally have much smaller growth rates than the SW modes, 
thereby better approximating the free modes' zero growth 
rates. Both the LW and SW modes of the asymmetric case 
exhibit shorter periods than those of the corresponding 
modes of the symmetric case. 

It is also noted that varying dr\/dz produces little 
change in frequency. Only a slight decrease in period is 
noticeable with an increase in 9r|/3z. The deeper boundary 
layer does not appreciably affect the wave period either, 
although a very small increase in frequency is detectable 
for the longer wavelengths of the symmetric LW mode and 
a small decrease for the asymmetric LW and SW modes, as 
well as for the symmetric SW mode. 

D. EIGENFUNCTIONS 

Figures 10 -and 11, which represent the free mode pro- 
files for = i and 2 km, respectively, are consistent 

with the solutions in Chang (1973b) . When CISK heating 
is included, the p and w^ profiles for both symmetry cases 
generally tilt nor theas t- sou thwes t between the northern 
boundary and the equator. However, the equatorward limit 
for the tilt in the p profile of the symmetric case depends 
upon wavelength (the longer the wavelength the farther away 
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from the equator is the limit). Such a northeast-southwest 
tilt in the eigenfunctions implies poleward transport of 
westerly momentum and agrees with observations summarized 
in Wallace (1971) . 

Comparison of Figures 12 and 16 (among others) shows 
that the tilt for the asymmetric LW mode decreases with 
longer waves. The SW mode generally shows a similar change. 

By increasing the vertical heating gradient at the top 
of the boundary layer, the tilt increases for both asymme- 
tric modes especially in the short waves (see Figures 12- 
15) . The only significant change in tilt due to an increase 
in the boundary-layer depth occurs for the LW mode, for 
which the tilt increases with decreasing wavelength (compare 
Figures 12 with 14 and 13 with 15) . 

Figures 12 and 16 show a decrease in the SW mode's 
equatorward maximum of w T amplitude and a general smoothing 
of the W T profile when the wavelength decreases. On the 
other hand, the LW mode shows a slight equatorward shift 
of the w^ maximum and a uniform amplitude increase with 
decreasing wavelength. The p profiles for both modes re- 
main nearly unchanged in these figures. 

An increase in 3r|/3z shifts the asymmetric LW mode's 

maximum in w and p more equatorward, particularly for the 
T 

shorter waves (see Figures 12 and 13) . The SW mode, mean- 
while,^ shows a slight equatorward shift only in its pole- 
ward maximum and a somewhat smoother profile of w T when 
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3r|/3 z increases. The p profile of the SW mode remains al- 
most unchanged with a larger 3r|/3z. 

Deepening the boundary layer also effects an equator- 
ward shift in the asymmetric LW mode’s p and w^ maxima 
(compare Figures 12 with 14 and 16 with 18). This result 
is similar to that found by Chang (1973b). From Figures 
16 and 18, it is apparent that for longer waves the SW 
mode’s equatorward maximum in w T becomes the primary maxi- 
mum when the bou nd ary - 1 ay e r depth increases. While both 
W T maxima of the SW mode shift slightly toward the equator 
with a deeper boundary layer, the p profile remains 
virtually unchanged (Figures 16 and 18) . Both modes indi- 
cate a larger amplitude when increases from 1 km to 

2 km, while the p amplitude again remains nearly constant 
(compare Figures 12 with 14 and 16 with 18). 

Essentially the same results as those just noted for 
the asymmetric case are found for the symmetric case. 

Only the significant differences are described below. 

From Figures 24 and 25 it is seen that the equatorward 
maximum of the symmetric SW mode shifts poleward for 
the longer wavelengths when 3r|/3z increases. In contrast, 
the asymmetric SW mode's equatorward w^ maximum remains 
stationary when 3p/3z is increased but the poleward maxi- 
mum moves equatorward. While an increase in the boundary- 
layer depth shifts the asymmetric LW mode's and p maxima 

more toward the equator, only the w maximum moves 
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equatorward for the symmetric LW mode (compare Figure 20 
with 22 and 24 with 26). 

The phase difference between the pressure and vertical 
motion fields generally remains 1/2 cycle for both symme- 
try cases under the conditions investigated. The only 
exception appears with the symmetric case. In the vicinity 
of the critical latitude and equatorward, a rapid phase 
shift in w„ occurs such that the two fields are nearly in 
phase at the equator. Whereas an out-of-phase relation- 
ship implies a downward energy transfer from the source 
region immediately above the boundary layer, in-phase im- 
plies the converse transport of energy, that is, from the 
boundary layer to the atmospheric layer immediately above. 
However, the relative amplitudes of both p and w are 
small near the equator, thereby minimizing the upward 
transport of energy within the ne ar - equa t or ia 1 region. 
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V. SUMMARY AND CONCLUSIONS 



A linear, tropical boundary-layer model is formulated 
on an equatorial beta-plane by utilizing the primitive 
equations of motion, a simplified thermodynamic energy 
equation and a CISK parameterization. Solutions in the 
form of complex wave frequencies are obtained numerically 
as eigenvalues to a second-order finite difference equa- 
tion . 

The primary conclusion reached is that there is no 
preferred sy nopt i c - s ca 1 e wavelength for CISK-forced tropi- 
cal waves. Although two modes for each symmetry case 
exist when the CISK forcing is incorporated., in no case 
does a noticeable peak in the growth rate curves 
materialize in the synoptic scale range. In addition, 
the e-folding times for the more realistic vertical heat- 
ing gradient (dr\/dz = IS ) are found to be of the order of 
5 days or slower, even for the deeper boundary- layer 
depths. This magnitude of growth rate is quite small 
compared to baroclinically unstable waves in the midlati- 
tudes. Consequently, even if one scale is slightly pre- 
ferred among the others, its relatively faster growth 
would be at such a slow rate that the competition for 
energy among the waves may not be strong enough to pro- 
duce a dominant wavelength. 
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The main conclusion of this study - that the growth 
rate is nearly uniform in the synoptic scale range - agrees 
with some earlier works (Chang, 1971, and Williams and 
Robertson, 1973), but differs from others (Yamasaki, 1969, 
1971; Hayashi, 1971;and Murakami, 1972). A similar inde- 
pendence of scale is noted for the so-called "spin-down" 
problem. By specifying 3ri/3z = -1 and substituting it into 
the forcing function on the right-hand side of .(15), a 
crude formulation of the spin-down problem is made. By so 
doing, the vertical motion profile for the free atmosphere 
is assumed to decrease linearly from w^ at the top of the 
boundary layer to zero at one scale height (i.e., z = H). 
The resultant damping rates are found by using the same 
procedure for computing the CISK growth rates and they are 
shown in Figure 28. The spin-down times are small as 
expected because of the small rotation rates near the 
equator. In addition the curves in Figure 28 show nearly 
uniform spin-down rates for the synoptic-scale waves. The 
CISK model equation, therefore, is mathematically analogous 
to the spin-down equation except that the forcing function 
has the opposite sign. Hence, it may not be surprising 
that the CISK growth rates are also independent of scale. 

Although there is no preferred wave scale in the 2000- 
10,000 km range, the LW mode of both symmetry cases never- 
theless exhibits larger growth rates in the synoptic scale. 
In contrast, the short wave end of the LW mode has decreas- 
ing growth rates with decreasing wavelength. This mode 
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with its growth rate characteristics seems to agree with 
the HB mode of Yamasaki (1969) and the findings of Murakami 
(1972). In addition, their results are similar to the LW 
mode in the following two ways: 1) maximum amplitude 
exists near the equator, and 2) wave periods of about 5 
days occur. On the other hand, the SW mode, especially 
for the shallower boundary layers, has its maximum ampli- 
tude more poleward and possesses larger periods. The 
SW mode also appears to prefer zero scale for growth. 

These three characteristics indicate that the SW mode may 
be identified as a hurricane (or typhoon) mode. 

From the results of this study, a growth rate depen- 
dence upon boundary- lay er depth appears to be two-fold. 
First, by increasing the depth of the boundary layer, 
growth rates increase for both the LW and SW modes of both 
symmetry cases. Secondly, the greatest change in the 
growth rates for both cases generally occurs between z 

= 1 km and z =2 km. As a result of these findings, it 
T 

is apparent that linear, tropical boundary- layer models 
must account for the variations in boundary- lay er depth 
in order to describe more accurately the growth rates of 
tropical lower tropospheric waves under CISK forcing. 

Within the tropical boundary layer, there is generally 
a transition zone between a quasi-Ekman boundary layer of 
the more poleward latitudes and a quasi-Stokes boundary 
layer of latitudes near the equator. The s low-rotation 
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Stokes boundary layer with its inefficient pumping does not 
contribute much to wave growth. The inclusion of the tem- 
poral acceleration terms in the bou ndar y- 1 ay er equations 
thus seems to have very little effect in CISK instability. 
In fact the results of this study generally agree with the 
qua s i -geos trophi c models of Chang (1971) and Williams and 
Robertson (1973). 

At the top of the boundary layer, the vertical motion, 
w , assumes a smoother profile across the transition zone 
(critical latitude) for the growing solution compared to 
the free solution in Chang (1973b). However, the large 
phase shift across the transition zone remains similar to 
that found for the free solution. 

Non-linear advective terms may be important in affect- 
ing the growth rates of synoptic-scale tropical waves. 

Such terms are neglected in this study and should be in- 
cluded in future research to determine their effect in a 
CISK model. Finally, a possible improvement to this linear 
model may be the 'inclusion of lateral shear of the mean 
flow in the equations for both the boundary layer and the 
free atmosphere. 
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Figure 1- Wave periods for free modes. 
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Figure 2. Growth rates for asymmetric case and = 1 km. 
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Figure 4. Growth rates for symmetric case and = 1 km. 
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Figure 5. Growth rates for symmetric case and z T = 2 km. 
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Figure 6. Wave periods for asymmetric case and z T - 
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Figure 7. Wave periods for asymmetric case and z^ - 2 km. 
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Figure 8. Wave periods for symmetric case and z = 1 km. 
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Figure 9. Wave periods for symmetric case and z ^ — 2 km. 
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Figure 10. Eigenfunctions for free mode of 
case (two wavelengths given), 
indicate pressure perturbation, 
lines, the vertical motion, w f 
the boundary layer. 
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Figure 11. Eigenfunctions for free mode of symmetric 
case (two wavelengths given) . Solid lines 
indicate pressure perturbation, p, and 
dashed lines the vertical motion, w^, , at the 
top of the boundary layer. 
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Figure 12. Eigenfunctions for 2000 km wavelength, 

asymmetric case, z = 1 km and 3n/3z = Is. 

Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 

w , at the top of the boundary layer. 
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Figure 13. Eigenfunctions for 2000 km wavelength, 

asymmetric case, z - 1 km and 9g/9z = 2S . 
Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 
w , at the top of the boundary layer. 
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Figure 14- Eigenfunctions for 2000 km wavelength, 

asymmetric case, 2 = 2 km and dr\/dz = IS . 

Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 

w , at the top of the boundary layer. 
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Figure 15. Eigenfunctions for 2000 km wavelength, 

asymmetric case, z = 2 km and 3f|/3z = 2S - 

Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 
w , at the top of the boundary layer. 
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Figure 16. Eigenfunctions for 4500 km wavelength, 

asymmetric case, z = ^ km an< ^ 3ri/8z = Is. 

Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 

w . at the top of the boundary layer. 
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Figure 17. Eigenfunctions for 4500 km wavelength, 

asymmetric case, z = 1 km and 8ri/3z = 2S . 

Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 
w T , at the top of the boundary layer. 
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Figure 18. Eigenfunctions for 4500 km wavelength, 

asymmetric case, z = 2 km and 8f|/3z = IS . 
Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 

w , at the top of the boundary layer. 
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Figure 19. Eigenfunctions for 4500 km wavelength, 

asymmetric case, z = ^ km an< ^ 3n/3z = 2S . 
Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 
w , at the top of the boundary layer. 
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Figure 20. Eigenfunctions for 2000 km wavelength, 

symmetric case, z = 1 km and 3n/9 z = IS. 
Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 
w , at the top of the boundary layer . 
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Figure 21. Eigenfunctions for 2000 km wavelength, 

symmetric case, z = 1 km and 3n/3z = 2S . 
Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 
w , at the top of the boundary layer. 
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Figure 22. Eigenfunctions for 2000 km wavelength, 

symmetric case, z = 2 km and 8n/9z = IS. 

Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 
w , at the top of the boundary layer. 
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Figure 23. Eigenfunctions for 2000 km wavelength, 

symmetric case, Z T = 2 km and dl)/dz = 2S . 
Solid lines indicate pressure per tr uba t i on , 
p, and dashed lines the vertical motion, 

w # at the top of the boundary layer. 
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Figure 24. Eigenfunctions for 4500 km wavelength, 

symmetric case, z - 1 km and dr]/3z = IS. 

Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 

w , at the top of the boundary layer. 
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Figure 25. Eigenfunctions for 4500 km wavelength, 

symmetric case, z = 1 ^m an< ^ 9n/9z = 2S . 

Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 

w , at the top of the boundary layer. 
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Figure 26. Eigenfunctions for 4500 km wavelength, 

symmetric case, z = 2 km and dr\/dz = IS . 

Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 

w , at the top of the boundary layer. 

T 



56 



LATITUDE 



W 



T 




Figure 27. Eigenfunctions for 4500 km wavelength, 

symmetric case, z rp = 2 km and dr\/dz = 2S. 
Solid lines indicate pressure perturbation, 
p, and dashed lines the vertical motion, 

w , at the top of the boundary layer. 
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Figure 28. Spin-down times for z ^ ~ 2 km. 
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